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We study the photon-triggered light and heavy meson production in both p+p and A+A collisions. 
We find that a parton energy loss approach that successlully describes inclusive hadron attenuation in 
nucleus-nucleus reactions at RHIC can simultaneously describe well the experimentally determined 
photon-triggered light hadron fragmentation functions. Using the same framework, we generalize our 
formalism to study photon-triggered heavy meson production. We find that the nuclear modification 
of photon-tagged heavy meson fragmentation functions in A+A collision is very different from that of 
the photon-tagged light hadron case. While photon-triggered light hadron fragmentation functions 
in A+A collisions are suppressed relative to p+p, photon-triggered heavy meson fragmentation 
functions can be either enhanced or suppressed, depending on the specific kinematic region. The 
■ anticipated smaller energy loss for 6-quarks manifests itself as a flatter photon-triggered B-meson 

fragmentation function compared to that for the D-meson case. We make detailed predictions 
for both RHIC and LHC energies. We conclude that a comprehensive comparative study of both 
photon-tagged light and heavy meson production can provide new insights in the details of the jet 
quenching mechanism. 
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High transverse momentum partons and their hadronic fragments are powerful and valuable probes of the high 
energy density matter created at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). 
These energetic partons are created in the early stage of the collisions and therefore provide access to the space- 
, time history of the transient hot and dense nuclear medium - the quark-gluon plasma (QGP) created in the heavy 
^\ • ion reactions. Specifically, as they propagate through the QGP, they interact with the medium and lose energy, a 
phenomenon known as "jet quenching" [l|-[8|. 

Experimentally, there has been tremendous progress in recent years [9[ in establis hing th e jet quenching phenomenon 
\Q | from different observables, such as the suppression of single hadron production |10Ml2j . dihadron correlations 1131 . 
[HI , 7-hadron correlations [lfl [l6[ , and ultimately the alteration of inclusive jet [l7j and dijet production [3 Il9j |. 
Theoretically, many perturbative QCD-based models of jet quenching have been developed and used to extract 
the medium properties [201 ] . especially the QGP parton number and energy densities. They are able to successfully 
describe most of the experimental data and the measured nuclear modification factor Raa of single hadron production 
in particular. 

To study jet quenching and parton energy loss in more detail and thus constrain the medium properties better, 
more sophisticated observables have been proposed to improve our understandin g of inelastic parton interactions in 



the QGP. One such example is 7-triggered hadron production and correlations |21|: one studies jet quenching by 
measuring the px distribution of charged hadrons in the opposite direction of a trigger direct photon. Since the 
photon does not interact with the dense medium, its energy approximately reflects that of the initial away-side parton 
before energy loss up to corrections at next-to- leading order in a s . One can therefore study the effective medium 
modification of the jet fragmentation function. Earlier studies on the 7-triggered light hadron correlations (22T - [2^ | 
seem to be roughly compatible with the experimental data. 

Although jet quenching has been successful in describing most of the experimental findings, there are still remaining 
puzzles. Resolving these puzzles will eventually uncover the exact underlying mechanism for the suppression of leading 
particles and jets. One well-known difficulty is related to the fact that the c and 6-quark parton-level energy loss 
in the QGP has not been sufficient in the past to explain the large suppression of non-photonic e + + e~ measured 
at RHIC [25], HH- These non-photonic electrons come from the semileptonic decays of D and B mesons. In the 
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framework of perturbative QCD, collisional dissociation of heavy-mesons in the QGP has been suggested as an 
additional suppression mechanism [27], [28j]. To distinguish with confidence between theoretical models one hopes 
to have a direct measurement of the heavy meson cross sections in both p+p and A+A collisions. Such direct 
measurements are finally becoming available from both RHIC and the LHC experiments. 

The main goal of this paper is to investigate new experimental observables, such as 7-triggered hadron production, 
and make simultaneous predictions for both light and heavy meson final states. Naturally, we will first focus on the 
standard parton energy loss mechanism in the QGP. More specifically, we present detailed studies of both 7-triggered 
away-side light hadron and heavy meson spectra in heavy ion collisions. Within the same energy loss formalism - the 
Gyulassy-Levai-Vitev (GLV) approach - we study how the QGP medium affects the 7-triggered light and heavy meson 
effective fragmentation functions. We find that the nuclear modification factor Iaa behaves very differently depending 
on the parent parton mass: Iaa is considerably suppressed for light hadrons {Iaa < 1), whereas it can be suppressed 
(Iaa < 1) or enhanced (Iaa > 1) for heavy mesons, depending on the kinematic region, due to the different shape 
of the heavy quark fragmentation functions. Furthermore, we find that the nuclear modification Iaa is flatter in the 
B-meson case compared to that of D-mesons due to the smaller energy loss for 6-quarks in the non-asymptotic (finite 
energy) case. Thus, a comparative study of 7-triggered light and heavy meson correlations can be a very useful probe 
of the physics that underlays the experimentally established jet quenching. 

The rest of our paper is organized as follows: in Sec. II we present the relevant formalism for both photon-tagged 
light hadron and heavy meson production. In Sec. Ill we present our phenomenological studies. We first compare our 
calculation to the experimental data on photon+light hadron correlations. Then, we make predictions for photon- 
triggered light and heavy meson production relevant to both RHIC and LHC experiments and discuss the differences 
in the observed nuclear modification. We conclude our paper in Sec. IV. 



II. PHOTON- TAGGED LIGHT AND HEAVY MESON PRODUCTION 



In this section we present the relevant formalism for the photon-triggered light and heavy meson production in 
nucleon-nucleon collisions. These formula will then be used in the phenomenological studies in the next section. 



A. Photon-tagged light hadron production 



Within the framework of the collinear perturbative QCD factorization approach, the lowest order (LO) differential 
cross section for back-to-back photon-light-hadron production can be obtained from the partonic processes: qq — > 717 
and qg — > jq and is given by [29j 
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where S is the squared center of mass energy of the hadronic collisions, and zt, x a , and Xb are given by 
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Here, y 7 and px (yh and Pt h ) are the rapidity and transverse momentum of the photon (away-side hadron), re- 
spectively. We denote by f a ,b/N(x a ,b) the distribution functions of partons a,b in the nucleon and Dh/ C (z) is the 

fragmentation function of parton c into hadron h. \M\ ab ^ are the squared matrix elements for ab — > 7c partonic 
processes, and are given by 
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where s, t, u are the partonic Mandelstam variables, e q is the fractional electric charge of the light quark, and N c = 3 
is the number of colors. 

Experimentally, one typically defines the so-called 7-triggered fragmentation function 



d aw( z t) = 
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for nucleon-nucleon collisions and a similar D"^ a (zt) per binary scattering for A+A collisions. The denominator 
defines the normalization and is given by 

= W dVHd PTh , .*f» , , (6) 
dy 7 dpT-, J dy 7 dyh dpr y dp Th 

which is just the cross section for direct photon production. To quantify the modification of 7-triggered fragmentation 
function in A+A collisions relative to that in nucleon-nucleon collisions due to the jet quenching, one introduces the 
nuclear modification factor, 

frh ( s _ d aa( z t) , 7 s 
d nn( z t) 

Note that, in spite of their suggestive name, 7-triggered fragmentation functions are derived from the observed away- 
side meson distribution and thus reflect all input in a pQCD calculation - the parton distributions, the hard scattering 
cross sections, the medium-induced radiative correction (or parton energy loss) and the parton decay probabilities - 
not only the fragmentation process itself. 



B. Photon-tagged heavy meson production 



Within the perturbative QCD factorization approach there have been different ways to calculate heavy quark 
production. One of them is called a fixed-flavor-number scheme (FFNS) [30l - [3^ | and is based on the assumption 
that the gluon and the three light quarks (u, d, s) are the only active partons. The heavy quark appears only in the 
final state and is produced in the hard scattering process of light partons. The heavy quark mass m is explicitly 
retained throughout the calculation. The other approach is called a variable-flavor-number scheme (VFNS) |33l - [35| 
and describes the heavy quark as a massless parton of density fg/j^ix, /1 2 ) in the nucleon, with the boundary condition 
/q/at(x, [i 2 ) = for [i < m. Thus, the heavy quark mass m is set to zero in the short-distant partonic cross section 1 . 

In this paper, we are particularly interested in studying the jet quenching effects for both light and heavy mesons. 
The different amount of energy loss during quark propagation in the medium created in heavy ion collisions is due 
to the mass difference between light and heavy quarks, at least in perturbative QCD. With this in mind, we want to 
keep explicitly the heavy quark mass in our calculation. In other words, we favor the FFNS scheme for the project at 
hand. Within this scheme, the photon-tagged heavy mesons are produced through the following partonic processes 
- (a) quark-antiquark annihilation: qq — > QQ^f; (b) gluon-gluon fusion: gg — > QQ'y. The sample Feynman diagrams 
are given in Fig. [TJ It is important to realize that in the FFNS scheme photon+heavy quark events are generated to 
LO by the 2 — > 3 processes, to be compared to the usual 2 — > 2 processes for the photon+light hadron events at LO. 






FIG. 1: Sample Feynman diagrams at leading order for direct photon plus heavy quark production: (a) for light quark-antiquark 
annihilation qq — > QQ"/, (b) and (c) for gluon-gluon fusion gg —¥ QQ'y. 

The differential cross section for photon-tagged heavy meson production takes the following form (32|: 
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1 Strictly speaking, this is the so-called zero-mass variable-flavor-number scheme (ZM-VFNS). There is also general-mass variable-flavor- 
number scheme (GM-VFNS) which combines the virtues of the FFNS and the ZM-VFNS, though more complicated. 
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where Dh/q{z) is the heavy quark Q to heavy meson H fragmentation function and <f> = <f>h — <t>i € [ 7r /2, 37r/2] is 
the azimuthal angle between the triggered photon and the associated heavy meson on the away side. We denote by 

totq — \ Pt + m2 the transverse mass and by pt q = Pt h I ' z the transverse momentum for the heavy quark. The 



(9) 



parton momentum fraction Xb is fixed by the kinematics and has the following form 
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The partonic matrix elements \M\ ab ^ > QQ 1 can be calculated in perturbative QCD. Note that the cross section 
does not produce any singularity since the heavy quarks have been taken to be massive explicitly. Therefore no 
regularization is needed despite the fact that we are act ually dealing with a 2 — » 3 process. The cross section for 
photoproduction of heavy quarks has been published in [36]. By using crossing symmetry, we can infer the spin- 
averaged matrix elements in our case. We list the results here for completeness. For the quark-antiquark annihilation 
channel q(px) + q{p 2 ) -t Qips) + Q(pa) + lips) we have 
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where eQ is the fractional electric charge of the heavy quark, and the coefficients Ax, A 2 , and A 3 are given by 
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Here, we have used the notation pij 
Qips) + QiPi) + lips) we have 
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The Rqed has to be summed over 12 permutations corresponding to the 6 permutations of the momenta px, P2, ~P5 
and two permutations of the momenta px and p2 ■ Rkf has to be summed over 4 permutations corresponding to the 
interchange of px, P2 and p 3 , pa- One such expression for Rqed is given by 
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and Rk f is given by 
Rkf 
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With the cross sections at hand, the experimentally accessible 7-triggered heavy meson fragmentation function is 
defined as 
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where the denominator is the differential cross section of photon and heavy meson production summed over all the 
heavy mesons 
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which is just part of the inclusive photon production cross section that has 7 + Q events. 

Likewise, we can define a 7-triggered heavy meson fragmentation function D^^(zt) in A+A collisions and the 
nuclear modification factor IJ^(zt) — D aa (zt) / DJJ^(zt) ■ In the next section we will study how these fragmentation 
functions behave in both p+p and A+A collisions and, most importantly, how they are modified in A+A collisions 
due to jet quenching. 



III. PHENOMENOLOGICAL STUDIES IN P+P AND A+A COLLISIONS 

In this section we will study the 7-triggered light hadron and heavy meson production in both p+p and A+A 
collisions. Final-state medium-induced radiative corrections are process dependent - they depend on the parameters 
of the QGP medium and, consequently, on the details of the heavy ion collision. However, in QCD they factorize from 
the hard scattering cross section [8j and enter the physical observables as a standard integral convolution. From the 
point of view of phenomenological applications, it is convenient to change the order of integration and include parton 
energy loss as an effective replacement of the fragmentation function _D^/ C (z) j2Z|: 

D h/c {z) => J de P(e) -^—D h /c {jZT^j • ( 20 ) 
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Here, P(e) is the probability distribution for a parton to lose a fraction of its energy e = ^ i due to multiple 
gluon emission. In our calculation P(e) is obtained using the GLV formalism for evaluating the medium-induced 
gluon bremsstrahlung [H, [3!| and the Poisson approximation. Note that it is not mandatory to employ Eq. ([20]) in jet 
quenching phenomenology. One can alternatively evaluate the attenuated partonic cross section prior to fragmentation 
and obtain identical results [28[ . 

We take into account parton energy loss in 7-hadron correlation in high energy heavy-ion collisions as in Eq. (|20j) for 
simplicity. We do not perform a new evaluation of the energy loss since we aim at a direct and consistent comparison 
between different final states in the same jet quenching formalism. Instead, we employ the differential medium-induced 
bremsstrahlung distributions and energy loss probabilities P(e) for light partons and heavy quarks obtained in [28ll40|. 
An illustration of these probabilities 

J! p{e)de = h J! ePi£)d6 = (^r)> e = £^r (21) 

is given in Fig. [2] Note that, in the probabilistic application of parton energy loss, larger (AE/E) corresponds 
to a distribution skewed toward large e and smaller (AE/E) corresponds to a distribution skewed toward small e. 
All results include an average over the jet production points distributed according to the binary collision density in 
an optical Glauber model. These energetic jets propagate through the Bjorken expanding medium that follows the 
number of participants density [28|, |40j. The left panel of Fig. [5] corresponds to central Au+Au collisions at RHIC 
and parton energy E = 10 GeV and the right panel of Fig. [5] corresponds to central Pb+Pb collisions at the LHC and 
parton energy E = 25 GeV. 

The difference between light quarks and gluons arises from the quadratic Casimir, or average squared color charge, 
of the parton. Note that the naive relation (AE) g = CA/CF{AE) q is approximately fulfilled only if (AE) q>g -C E. 
In the realistic case, (AE) g < C^/Cf(A£) ? and the deviation form the naive scalin g ca n be significant. On the 
other hand, the difference between light and heavy quarks is related to the quark mass |25l. l2ol |39| . In the very high 
energy limit, for coherent medium-induced bremsstrahlung, the mass dependence of parton energy loss disappears. 
The calculations that we present in this paper are not in this asymptotic limit. Furthermore, for partons of effective 
mass gT ~ m c , where g is the coupling constant and T is the temperature of the quark- gluon plasma created in heavy 
ion collisions, there is very little difference between the light quarks and the charm quarks. This is clearly seen in 
both panels of Fig. [2] and well understood. For the much heavier bottom quarks m& 3> m c ~ gT, however, the energy 
loss is noticeably smaller and P(e) is skewed toward smaller values of e. 




FIG. 2: The probability distribution P(e) for partons to lose a fraction e of their energy is shown for gluons (dot-dashed), light 
quarks (dashed), charm quarks (solid) and bottom quarks (dotted). Left panel: central Au+Au collisions at RHIC and parton 
energy Ej — 10 GeV. Right panel: central Pb+Pb collisions at the LHC and parton energy Ej = 25 GeV. 



A. Comparison to experimental data: 7-triggered light hadron production 

In this subsection we will compare our calculations for 7-triggered light hadron production to the experimental data 
at RHIC. We will use CTEQ6 parton distribution function 4l| and deFlorian-Sassot-Stratmann (DSS) light hadron 
fragmentation function [42j. 

In Fig. [3] (left and top right panels) we show our results compared to both PHENIX [l5[ and STAR data [l|| ■ For 
p+p collisions we use the definition given in Eq. ([5])- As seen from these figures, the perturbative QCD calculation gives 
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a very good descriptions of the experimental data. For Au+Au collisions we use the same energy loss distributions 
that were used to describe single hadron suppression. As can be seen from Fig. |31 the calculated D a h a {zt) are in 
good agreement with the experimental data with the possible exception of the large zt — > 1 region in the STAR 
data. This part of phase space is also inherently difficult to measure accurately. Such agreement gives an independent 
confirmation of the jet quenching mechanism in light hadron production. 
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FIG. 3: 7-triggered fragmentation functions D(zt) (left) and (top right) are plotted as a function of zt at ^/Snn = 200 GeV. 
Left: top for p+p collisions, bottom for Au+Au collisions. Photon and hadron rapidities are integrated over [—0.35, 0.35], while 
the trigger particle (the photon) momentum has been integrated over [5, 7], [7, 9], [9, 12], and [12, 15] GeV from top to bottom. 
Data is from PHENIX [15|]. Right: top for D(zt) - solid for p+p collision, and dashed for Au+Au collision; bottom for nuclear 
modification factor I a A- The photon and hadron rapidities are integrated over [—1,1], and the photon momentum is integrated 
over [8, 16] GeV. Data is from STAR 
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FIG. 4: Nuclear modification factor Iaa(zt) is plotted as a function of zt at ^/Snn = 200 GeV. The numbers on the right 
top corner in each plots are the range of the photon momentum, while the photon and hadron rapidities are integrated over 
[-0.35,0.35]. The data is from PHENIX [lj|. 



It is worthwhile to point out that the 7-triggered light hadron fragmentation function D aa (zt) in Au+Au collisions 
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is suppresses relative to that in p+p collisions. This is understandable since at LO D lh (zr) oc £>/j/ c ( z t) according 
to Eq. (fTJ) , where D h / C (z) is a fast falling function of z. Since jet quenching can effectively be thought of as a 
mechanism that probes slightly higher z as in -Dh/ C ( jzr) according to Eq. (|20|) . sampling higher z leads to a smaller 
fragmentation function. Thus, the 7-triggered fragmentation function in A+A collisions will be suppressed compared 
to p+p collisions. We will find this behavior can be altered in the 7-triggered heavy meson production as we will 
demonstrate below. 

In Fig. [3] (bottom right panel) and Fig. |4j we compare our calculations of I AA to the experimental data. Being a 
ratio of the 7-triggered fragmentation functions in A+A and p+p collisions, the nuclear modification factor I aa (zt) — 
D aa (zt) I Dpp ( z t) has much larger error bars. Our results are consistent with the experimental data (within the error 
bars). 

B. Predictions for 7-triggered light and heavy meson production 

In this subsection we will make predictions for the 7-triggered light and heavy meson production at both RHIC 
and LHC energies and will compare the different features of the in-medium modification depending on the parton 
mass. We use D and J3-meson fragmentation functions as calculated in [28| and choose m — 1.3 (4.5) GeV for the 
charm (bottom) quarks. For 7-triggered heavy meson production we have integrated over the relative azimuthal angle 
4> = 4>h ~ 4>i m the region: \(f> — tt\ < n/5 in Eq. (jH)), following the experimentally applied cuts [TBI flo|. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

FIG. 5: Top panels: predictions for 7-triggered fragmentation functions D(zt) where the solid lines are for p+p collisions and 
the dashed lines are for central Au+Au collisions. Bottom panels: predictions for the nuclear modification factor Iaa(zt) 
where the solid lines are for B-meson, the dashed lines are for D-meson, and the dash-dotted lines are for charged hadrons. 
The left plot has kinematics similar to STAR, while the right plot has kinematics similar to PHENIX. 

In the upper panels of Fig.[5]we plot D(zt) as a function of zt at V Snn — 200 GeV for both p+p and central Au+Au 
collisions. Left panels reflect STAR kinematics and right panels reflect PHENIX kinematics. As seen clearly in the 
plots, D(zt) for light hadrons and heavy mesons are very different. For light hadrons, the 7-triggered fragmentation 
function D i}i (zt) drops very fast as zt increases, consistent with the behavior of light hadron fragmentation function 
D h / C (z). On the other hand, the 7-triggered heavy meson fragmentation function D lH (zt) is a relatively flat function 
of zt- Even though in the FFNS scheme D iH {zt) does not directly correspond to the heavy quark to heavy meson 
decay probability D h /q{z) itself, it docs reflect to some extent the major feature of the heavy meson fragmentation 
function: it grows at lower z and decreases at higher z. This difference in the fragmentation functions will lead to 
distinctive difference in the nuclear modification factor I AA between light hadron and heavy meson. 

In the bottom panels of Fig. [S] we plot Iaa (zt) as a function of zt at RHIC for 7-triggered light charged hadrons 
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FIG. 6: Top panels: predictions for the 7-triggered fragmentation functions D(zt), where the solid lines are for p+p collisions 
and the dashed lines are for central Pb+Pb collisions at ^/Snn = 2.76 TeV. Bottom panels: predictions for the nuclear 
modification factor Iaa(zt), where the solid lines are for B-meson, the dashed lines are for D- meson, and the dash-dotted 
lines are for charged hadrons. We have integrated the photon and hadron rapidities over [—2, 2]. For the left plot, the photon 
momentum is integrated over [10, 20] GeV, while for the right plot, it has been integrated over [25, 50] GeV. 

(dot-dashed), D- mesons (dashed) and i3-mesons (solid), respectively. They exhibit very different behavior. For light 
hadrons, one expects I\ A < 1 due to jet quenching, as explained above. The magnitude of this suppression arrises 
mainly from the steepness of the light hadron fragmentation function D^iJ^z). On the other hand, according to 

Eq. ([5]). for the 7-triggered heavy meson case IZ A really depends on the whole z-integration of the heavy meson 
fragmentation function Djj/q(z). If the integral is dominated by the small- z region, where Djj/q{z) grows with 
increasing z, jet quenching in Eq. (|20p means sampling relatively larger- z and thus larger D h /q(z). Consequently, 

one will then have > 1. On the other hand, if the integral is dominated by the large-z region, where D h /q(z) 

decreases with increasing z, sampling relatively higher-z means smaller Dh/q{z). One thus has I AA < 1. 

One keeps in mind that we have three-particle final state, thus zt = Pt h /pt~, is not the same as the momentum 
fraction z in heavy meson decay probability D h /q{z) according to Eq. ©. Nevertheless, we find that the average (z) 

in the collision does increase as zt increases. Thus, at high zt where the z-integral in I AA is dominated by the large 
z region, one should expect that both B and _D-mesons are suppressed due to jet quenching. The magnitude of the 
suppression should follow I AA > I AA > I AA if the heavy quark loses less energy than the light quark, as predicted 
by perturbative QCD calculations. On the other hand, in the low zt region, according to our calculation, we find 
that I AA > 1 for _D-meson, consistent with our naive expectation, that is the low z region dominates the z integral in 
Eq. ®. 

However, for 7-triggered S-mesons t^ A B A < 1 for the whole zt region for the kinematics we have chosen at RHIC. 
This is due to the fact that our £>-meson fragmentation function is even harder than the D-meson fragmentation 
function. It drops very fast at high z while increases only slowly at low z. Thus the nuclear modification from 
high z (suppression) wins over that from low z (enhancement) for the kinematic region we have chosen. Combining 
the analysis for both the low and the high zt ends, one immediately finds that the nuclear modification I AA for 
7-triggered B-meson fragmentation function is natter than that for D-meson case. Another important reason for this 
flatter behavior and the smaller size of the nuclear modification is the smaller energy loss of 6-quark in comparison 
to c-quark. Thus, the shape of the nuclear modification for the 7-triggered fragmentation functions can cary valuable 
information about the properties of medium-induced gluon bremsstrahlung. 

In Fig. [5] we give predictions for the 7-triggered light and heavy meson production in both p+p and central Pb+Pb 
collisions at y/SjyN — 2.76 TeV at the LHC. We integrate over both the photon and hadron rapidities from -2 to 
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2. In the left (right) panel, the trigger photon momentum is integrated from 10 < pr < 20 GeV (25 < pr < 50 
GeV). The uncertainty band comes from the fact that we have included a ~ 25% uncertainty in the magnitude of the 
energy loss in our jet quenching calculation. The behavior of Iaa follows from similar considerations for the relevant 
kinematics. It is also interesting to notice that for low both TJ^ and I'aa can be larger than unity in the chosen 
kinematic region. 

It will be illuminating and important to measure such correlations for both light and heavy meson production in 
p+p and A+A collisions at RHIC and at the LHC. 



IV. CONCLUSIONS 



We studied photon-triggered light hadron and heavy meson production in both p+p and A+A collisions. We 
found that the energy loss approach that was successful in describing single hadron production in A+A reactions at 
RHIC could simultaneously describe the STAR and PHENIX experimentally extracted photon-triggered light hadron 
fragmentation functions. Using the same theoretical framework, we generalized our formalism to study photon- 
triggered heavy meson production. To take into account the heavy quark mass explicitly, we followed the so-called 
fixed-flavor-number scheme and derived the differential cross section for photon+away-side heavy meson. We found 
that the nuclear modification of photon-tagged heavy meson fragmentation functions in A+A collision is very different 
from that of the photon-tagged light hadron fragmentation functions. This variance was determined to arise from the 
different shape of the decay probabilities for light partons into light hadrons and heavy quarks into heavy mesons, 
respectively. Comparing D and B-mesons, we predicted that the nuclear modification factor Iaa would be flatter 
for 7+£? production than the one for the 7+D case. This is directly related to the different amount of energy lost 
by heavy quarks with different mass. Thus, the different shape of Iaa in "f+h, 7+-D, and 7+-B productions can be 
a sensitive and quantitative probe of the strength of medium-induced parton energy loss. Finally, we made detailed 
predictions for both photon-triggered light and heavy meson at RHIC and at the LHC. We conclude by emphasizing 
once again that a comprehensive study of these new final-state channels will provide fresh insight into the details of 
the jet quenching mechanism. 
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